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Abstract 

Using about 3.9 million hadronic Z decays from e+e~ collisions recorded by the 
OPAL detector at LEP at centre-of-mass energies a/s ~ Mz, the branching ratio 

for the decay D^T r^Vr has been measured to be 

BR(D; ^ T-y^) = (7.0 ± 2.1(stat) ± 2.0(syst)) %. 
This result can be used to derive the decay constant of the D~ meson: 
fn^ = (286 ± 44(stat) ± 41(syst)) MeV. 
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1 Introduction 



The branching ratio of the purely leptonic Dg — > £ z/^ decayQ can be calculated [|l] 
using 




BR(D; ^ rv,) = ^rn^M (^1 - £f ) \Vj'T^il^, (1) 

where is the mass and the lifetime of the D~ meson, f^^ the D~ decay 
constant and Vcs the corresponding CKM matrix element. Gp denotes the Fermi 
coupling constant and the mass of the lepton. 

Several models for the calculation of the decay constant {-q^ exist: potential 
models predict fn^ in the range from 129 MeV to 356 MeV [ij, QCD sum rule 
models predict fo^ = (231 ± 24) MeV and lattice QCD calculations predict 
fc, = (240li) MeV l- 

The extraction of CKM matrix elements from B'^ — B oscillation measure- 
ments relies on these theoretical models for calculation of the decay constant for 
B mesons, /b, since a measurement of /b from B~ — > I'vi decays is currently not 
feasible. It is therefore important to measure fn^ to test the theoretical models 
used in the /b calculation. Measurements of fn^ in leptonic D~ decays have been 
pubhshed by WA75 [|, BBS §, E653 §, L3 0, CLEO g, and BEATRICE 0. 
The measured values lie between 190 MeV and 430 MeV. The current world 
average is 280 ± 48 MeV ^ . 



In this paper, we present a measurement of BR(Dg r z/^) and fn^ based 
on reconstruction of the decay sequence 

e+e" ^ Z ^ cc ^ D*-X 

U 7 Dg- 

^-lJllJr (i = e,fi). (^2) 

Only D~ — > T~Ur events from Z — > cc decays are considered, since a measure- 
ment of BR(Dj7 T'Ur) in Z — > bb events is systematically limited by the large 
uncertainty on the production rate of D~ mesons in Z — bb events. 

Hadronic r decays are difficult to distinguish from background and therefore 
only r decays into electrons or muons are used. Since the D~ mass cannot be 
reconstructed from a single particle in the final state, a neural network is trained 
on a preselected sample of hadronic Z events with one identified electron or muon, 
requiring the kinematics to be consistent with D~ — >• r~i/^ — > ^"z/^z/^z/^ decays. 

In the last step of the analysis D*^ 7Dg" decays are reconstructed in this 
Dj7 T^v-r enhanced sample by forming the invariant mass of the photon and 
the D~ candidate. This reduces the dependence on the Monte Carlo simulation 
of the background and increases the purity of the D~ sample. 



^Charge conjugate decays are implied throughout the paper. 
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Since the D^T — > t^ue decay is helicity suppressed, the r channel has the largest 
branching ratio of all leptonic channels. Eq. |I] predicts the branching ratio into 
electrons to be negligible, BR(D~ — »• e~z7e)/BR(D~ — >• T~i)r) < 10~^, due to the 
factor m^, whereas the branching ratio into muons, BR(Djr fJ'^i^n), is expected 
to be sizeable, BR(D~ fi~h'^)/BR{D^ — > r^Ur) = 0.103. Therefore the decay 
Dj7 —>■ fJ^'i^fi is included in the signal definition and the final result is corrected 
for this contribution. 



2 Detector, data sample and event preselection 



The OPAL detector is described in detail elsewhere |[TT|. Tracking of charged 
particles is performed by a central detector, consisting of a silicon microvertex 
detector, a vertex chamber, a jet chamber and z-chambers0. The central detector 
is inside a solenoid, which provides a uniform axial magnetic field of 0.435 T. The 
silicon microvertex detector consists of two layers of silicon strip detectors; for 
most of the data used in this paper, the inner layer covered a polar angle range of 
I cos 6*1 < 0.83 and the outer layer covers | cos 6*1 < 0.77, with an extended coverage 
for the data taken after the year 1996. The vertex chamber is a precision drift 
chamber which covers the range | cos^| < 0.95. The jet chamber is a large-volume 
drift chamber, 4.0 m long and 3.7 m in diameter, providing both tracking and 
ionization energy loss (dE/dx) information. The 2r-chambers provide a precise 
measurement of the ^-coordinate of tracks as they leave the jet chamber in the 
range | cos6\ < 0.72. 

The coil is surrounded by a time-of-flight counter array and a barrel lead-glass 
electromagnetic calorimeter with a presampler. Including the endcap electro- 
magnetic calorimeters, the lead-glass blocks cover the range | cos^| < 0.98. The 
magnet return yoke is instrumented with streamer tubes and serves as a hadron 
calorimeter. Outside the hadron calorimeter are muon chambers, which cover 
93% of the full sohd angle. 

For Monte Carlo studies, event samples have been generated using JETSET 
7.4 [1^ for multihadronic Z events and KORALZ 4.0 [jl3| for r pair events. 



Special signal samples have also been generated using JETSET. These consist of 
Z ^ cc — > D~X events with the decay sequence given in Eq. 2. The D^T —>■ t~i)j. 
signal is normalised using /(c D^r) = 0.121 [1^ and BR(Dg" t^J^t) = 7.0%. 



Tau polarisation effects are handled by the r decay library TAUOLA 2.4 

The data sample used in this analysis consists of about 3.5 million Z de- 
cays recorded during the period 1991-1995 and an additional 0.4 million Z events 
recorded for detector calibration purposes in 1996-2000. Events are only used 
if the silicon microvertex detector and the other main detector components rel- 



right handed coordinate system is used, with positive z along the e~ beam direction and 
X pointing towards the centre of the LEP ring. The polar and azimuthal angles are denoted by 
and 0, and the origin is taken to be the centre of the detector. 
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evant for the analysis were fully operational. Hadronic Z decays are selected 
based on the number of reconstructed tracks and the energy deposited in the 



calorimeter [jT6|. To ensure that the event is well contained within the accep- 
tance of the central detector, the polar angle of the thrust axis is required to 
satisfy | cos6't| < 0.8. 

Signal events are characterised by the presence of an electron or a muon from 
T decays and large missing energy. Electrons and muons are identified using 
neural networks [0, [1^ which are trained to identify leptons with a momentum 
greater than 2 GeV. Only events with exactly one identified electron or muon are 
selected. Electrons from photon conversions are rejected using a neural network 



conversion finder [19 



Each event is divided into two hemispheres by the plane perpendicular to the 
thrust axis. The hemisphere with less visible energy is required to contain the 
lepton. This hemisphere is selected to search for D~ T~Ur decays. To further 
enrich the sample in cc events, a loose anti-b tag ||20| is applied in the hemisphere 
opposite to the D^. Finally, at least 9 tracks are required in an event to reduce 
the background from Z ^ r+r^ events while keeping more than 97% of the signal 
events at this stage of the selection. 



3 Reconstruction technique 

A matching algorithm is used to avoid double-counting of particle momenta 
in the calorimeters and in the tracking detectors. The output of the matching 
algorithm - referred to as particles - are tracks and calorimeter clusters. 

If the missing energy in the event is only due to the neutrinos produced in 
the D~ decay, the energy and momentum of the D~ are exactly given by 

^D. = - E (3) 

j^lepton 

^D. = V~s- (4) 

j^lepton 

where ^/s is the e"'"e~ centre-of-mass energy. The summation is performed over 
all particles in the event except the lepton. The resulting mean reconstructed 
energy of the D~ is 27 GeV which is slightly larger than the true mean energy of 
26 GeV. 

Due to detector acceptance and resolution effects this method yields an energy 
resolution of 6.5 GeV and an angular resolution of 52 mrad where the resolution is 
defined as the sigma of a single Gaussian fitted to the distribution. To further im- 
prove the energy resolution, a kinematic fit is applied in which the energy and the 
absolute momentum of all particles (except the lepton) are varied independently 
from each other (i.e. varying their mass) using the constraint 
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El-Pl=Mj,^. (5) 

The values calculated from the deviations from the experimentally mea- 
sured values 

f ZT'fit 77'meas\2 

are minimized. This procedure yields an energy resolution of about 3.0 GeV. 
About 2% of the events are rejected because the kinematic fit does not converge. 
The efficiency e(Dj7 t^J^t) to reconstruct the Z — > cc ^ Y)*~X, D*^ — > 70^7 
'~fT~i?j. signal in the —>■ ^^v^Vt channel or in the e~z7gZ/^ channel at this 

stage of the analysis is about 30%. 



Selection of — > t u-r candidates 



s 



In the next part of the analysis a — T^Vr enriched sample is selected using 
neural networks. About 52% of the selected events used as input to the neural 
networks are expected to be Z ^ bb events, about 36% Z ^ cc events and the 
remaining 12% Z boson decays into light quarks (uds). The signal contribution is 
expected to be of the order 1%. For each channel (electron or muon) two neural 
networks are trained: one to separate signal from Z — > cc background events and 
one to distinguish between signal and Z — bb background. 

The light-quark background is not used in the training. Since D^T tv^ 
decays from b decays are not considered signal, they are included in the bb 
background. They constitute about 0.8% of the bb background events used as 
input to the neural network. 

The following variables are used in all four networks: 

• The reconstructed energy i?Ds of the D~ obtained from the kinematic fit 
(Fig. |l]a); the reconstruction method used for the energy and the momentum 
of the D~ is only valid for purely leptonic decays. Semileptonic background 
decays are expected to have a lower reconstructed D~ energy E-^^. 

• The lepton energy -Eiep; leptons in light-quark background events have on 
average lower energy than in signal events whereas leptons in background 
events from b — > £ decays have on average higher energy than in signal 
events due to the hard fragmentation of the b hadron. 

• The output of two additional neural networks trained to find b — ^ £ (Fig- lUb) 



and b ^ c — i> £ decays leptons originating from signal events have 
properties more similar to leptons from b — £ decays than from b — c — £ 
decays. 
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• The visible invariant mass determined from the tracks and clusters and the 
energy sum in the electromagnetic calorimeter (ECAL), both calculated in 
the D~ hemisphere; on average they are lower for signal events due to the 
energy carried by the neutrinos. 

The choice of input variables is optimized separately for each net. The follow- 
ing two variables are only used in the two neural networks rejecting bb background 
and in the neural network separating the muon channel from the cc background: 

• The momentum piep.Os of the lepton in the D~ rest frame for signal events; 
it is limited by the mass of the D~ meson to be piep,Ds < The 
Picp.Ds distribution is smeared by the experimental resolution. For leptons 
not originating from D~ decays this restriction does not exist, leading to a 
tail at higher picp,Ds- 

• The angle between the direction of the reconstructed D~ and the direction 
of the jet containing the D~; in signal events this angle is on average larger 
than in Z — s> qq background events. The jets are reconstructed by combining 
all particles - including the lepton - using a cone algorithm |^3[. The jet 
direction is then calculated excluding the lepton. 

Variables sensitive to the flavour of the event are used in the neural networks 
separating signal from cc background: 

• The highest momentum Pmax of any particle with a charge opposite to that 
of the lepton in the D~ hemisphere (Fig. 0c); in signal events this particle 
should originate from the fragmentation of the c quark which produced the 
D~ meson. On average, it is therefore expected to have less momentum than 
the highest momentum charged particle in cc and light-quark background 
events. 

• The angle aiep,Ds between the direction of the lepton in the D~ rest frame 
and the direction of the reconstructed D~ in the lab frame (Fig. |l|d); for cc 
background events the aiep,Ds distribution is broad and it is peaked around 
7r/2 while in signal events aiep,Ds is closer to vr. 

In the neural networks which reject bb background the following variables are 
used: 

• The number of tracks and the number of clusters in the D~ hemisphere; 
the number of clusters and the b likelihood in the opposite hemisphere. 

Four selected variables which display good signal versus background separa- 
tion are shown in Fig. |l|. Data and background Monte Carlo distributions are in 
good agreement. 
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The output distributions of the neural networks for the muon channel are 
shown in Fig. The output distributions for the electron channel are similar. In 
Fig. ^, two-dimensional distributions of the outputs of the neural networks are 
shown. A cut at 0.85 on all outputs is applied to select a D~ — > T~Ur enriched 
sample with a signal efficiency of about e(D^ t^i^t) = 9%. 

The discrepancies between data and Monte Carlo are largest for the neural 
networks used to reject Z ^ cc background below 0.2, far away from the cut. 
Discrepancies between the neural network input distributions of the data and the 
Monte Carlo simulation are treated as systematic uncertainties, as discussed in 
Section |^. 



5 Dg reconstruction 

In the DjT t^Vt enriched sample, events with photon candidates in the D~ 
hemisphere are used to reconstruct the decay D*~ — > 7D~. The D~ signal can 
then be observed as a peak in the D*" mass region of the 7D~ invariant mass 
distribution. 

Since the neural networks have been trained to find D^T t^Vj. events and 
not specifically D*" tDjT events, additional cuts on the Dg" t^v^. enriched 
sample are required to reduce the background in the 7D~ invariant mass distri- 
bution: 



The b likelihood as given by the b tagging algorithm [20| in the hemi 



sphere has to be less than 0.5 to further suppress bb background. 

• Using energy and momentum conservation, the missing energy in the hemi- 
sphere is reconstructed from the visible energy -E^fg™' in the hemisphere, the 
invariant mass of all particles in the hemisphere, Mhemi? and in the opposite 
hemisphere, Mopp, via the relation: 

TTihemi V-^ , ^^-'hemi ^"^opp T-ihemi frrs 
^miss = ^ + ^vis • (7) 

The missing energy -^^1^' has to be larger than 15 GeV. 

These two cuts reduce the efficiency for the signal to 8%. 

The photon is found using information from the electromagnetic calorimeter 
as described in [^. This method assigns a weight to each photon candidate 



corresponding to the probability for it to stem from a real photon. To accept 
a photon candidate, this weight has to be larger than 0.6. Only events with 
exactly one such photon candidate in the D~ hemisphere are accepted. The 
distribution of the photon energy after all previously defined selection cuts is 
shown in Fig. ^. Data and Monte Carlo simulation are in reasonable agreement. 
Finally, is required to be greater than 2.3 GeV. For smaller photon energies the 
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shapes of the invariant mass distributions of the photon and the D~ candidate 
become similar for background and signal. The energy resolution for photons 
with > 2.3 GeV as determined by the Monte Carlo is about 300 MeV and the 
angular resolution is about 5 mrad. 



6 Results 

The distribution of the invariant mass m(7D~) of the photon and the D~ candi- 
date for the events satisfying all selection criteria is shown in Fig. ^ In the signal 
region, m(7D~) < 2.36 GeV, there are 24.5 ±2.8 background events predicted by 
the Monte Carlo. The number of background events is determined by requiring 
the expected number of Monte Carlo events to be identical to the number of data 
events after the lepton identification cuts described in Section 2. 

The most important branching ratios have been adjusted in the Monte Carlo 
using the values in . The uncertainty on the background is due to the limited 
number of Monte Carlo events. This number is subtracted from the data which 
yields 22.5 ± 6.9 signal events. The uncertainty is the statistical uncertainty of 
the data. 

The efficiency to reconstruct Z ^ cc ^ D^^X, D*^ 7T^z^r events 

in the r~ ^~i7^Vt- channel is 0.9% and in the erVeVr channel 0.7%. If 

the D~ decays directly into muons via D~ — > ^~v^l the efficiency is 0.6%. 

The branching ratio BR(Djr t^Vt) is extracted using 

br(d; ^ T-v,) ^^"'^'^ 



2iVz ■ i?c ■ /(c ^ D-) ■ Py(D^ Ds) ■ BR(Dr ^ 7Ds 
1 



X- 
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BR(D7^T-Pr) 

where Ai'cand is the number of background-subtracted candidates in the signal 



region, Nz the number of Z decays, -Rc = 0.1671 ±0.0048 [jT4| the partial width of 
the Z decaying into a pair of charm quarks, /(c D^r) = 0.121 ± 0.025 [|l^ the 
production rate of D~ mesons in charm jets, e(Dj7 t^Ut) the efficiency for the 
signal and e(D^ A* '^m) efficiency for reconstructing D*~ — > 7D~ — > 7yU~z/^ 
decays. As discussed in Section |I|, we use BR(Dg" yU^z/^)/BR(Dg" t^J^t) = 
0.103. 

P^/(D*,Ds) is the ratio of cs mesons produced in a vector state (D*) with 
respect to the sum of the pseudoscalar (Dg) and vector states. For non-strange 
D mesons, Pi/(D*,D) has been measured by ALEPH DELPHI ^ and 



OPAL [|8l. The averaged value is Pl/(D^D) = 0.61 ± 0.03 To extrapolate 
this ratio to Dg mesons, the effect of the decays of L = 1 D** resonances and quark 
mass effects need to be taken into account. D** resonances contribute only in the 
case of non-strange mesons. This effect was estimated by OPAL to be smaller 
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than the experimental uncertainty |^ and is therefore neglected. Applying the 
correction factor for quark mass effects from ^9] yields Pi/(D*, Dg) = 0.64 ± 0.05 
where the full size of the correction is included in the uncertainty. This value is 
consistent with the ALEPH measurement of Pi/(D*, D^) = 0.60 ± 0.19 [|| . 

Using Py(D*, Dg) = 0.64±0.05 and the values given in Table |I]for the branch- 
ing ratios yields 

BR(D7 ^ T-Ur) = (7.0 ± 2.1(stat)) %. (9) 



7 Systematic Uncertainties 

Systematic uncertainties arise from the uncertainties in the branching ratios, the 
Monte Carlo modelling, selection efficiencies and the detector resolution. The 
resulting systematic uncertainties are summarised in Table ^ and described in 
more detail below. 

External sources: The external values used in the calculation of the branching 
ratios are each varied within their uncertainties. 

Monte Carlo statistics The uncertainty on the background rate and on the 
efficiencies e(D^ t^t^t) and e(Dj7 l^^^n) due to the limited number of 
Monte Carlo events is counted as systematic uncertainty. 

Background: To account for uncertainties in the determination of the back- 
ground rate, the number of background events is also calculated using the 
sideband of the m{^'D^) distribution, defined by m{;~^Y)^) > 2.4 GeV. The 
difference between the standard analysis and the predicted background rate 
using the sideband is 0.1 background events in the signal region. This dif- 
ference is taken as systematic uncertainty. 

Background composition: About 55% of the background is due to combina- 
tions of Dj7 candidates with photons which do not originate from the same 
decay. The remaining background mainly consists of D*~ 'jD~ and 
D*° — > 7D° decays. No D~ — > rur decays from b decays are expected 
in the final sample. The most important branching ratios have been var- 
ied within their uncertainties using BR(D*^ 7D ) = 0.016 ± 0.004, 
BR(D*° ^ 7D°) = 0.381 ± 0.029, BR(D° e-X+) = 0.0675 ± 0.0029 and 
BR(D° l^~X+) = 0.068 ± 0.008 The corresponding change of the 



result is taken as systematic uncertainty. 

Fragmentation: To determine the effect of uncertainties in the Monte Carlo 
description of the fragmentation of b and c quarks, the distribution of the 
scaled hadron energy, xe = ^E^/ ^/s, is reweighted within the experimen- 
tal uncertainties for b quarks, (xe) = 0.702 ± 0.008, and for c quarks, 
(xe) = 0.484 ± 0.008 The largest of the variations observed using the 
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fragmentation functions of Peterson et al., Collins and Spiller, and Kartvel- 
ishvili et al. |^ is taken as systematic uncertainty. 



Lepton spectrum: The exact shape of the lepton momentum spectrum for 
background events is not known. Therefore b— b— s>c^£ and c — £ 
decays are reweighted to reproduce the lepton momentum spectrum in the 



rest frame of the b or c hadron as predicted by the ACCMM |31], the 



ISGW |^2| and the ISGW** models. The same parameters as in [|T4 



are used. The largest difference between the results obtained using the 
different models is taken as systematic uncertainty. 

Tracking resolution: To take into account uncertainties in the modelling of 
the tracking resolution by the Monte Carlo, the reconstructed Monte Carlo 
track parameters are smeared by ±10% and the analysis is repeated. 
The largest difference between the results is taken as an estimate for this 
source of systematic uncertainty. 

Photon energy: The analysis is redone varying the photon energy scale in the 
Monte Carlo simulation by ±2% and the difference between the re- 



sults is taken as systematic uncertainty. Furthermore it was checked that 
the result obtained with the low purity sample in the range 1 GeV < 
< 2.3 GeV is statistically consistent with the result obtained for > 
2.3 GeV. 

Lepton identification efficiency: The electron identification efficiency has been 



studied in |Q and has been found to be modelled correctly within 4%. The 
muon identification efficiency has been studied in |l^, giving an uncertainty 
of 5%. 

Neural networks: Each of the input distributions has been compared between 
data and Monte Carlo simulation. The simulated distributions are reweighted 
for each input variable in turn to agree with the corresponding data dis- 
tributions, and the analysis is repeated with the weighted events. The 
resulting differences in the measured branching ratio are added in quadra- 
ture to obtain an estimate of the systematic uncertainty due to modelling of 
the input variables. This includes the uncertainty due to the Monte Carlo 
modelling of the missing energy -Eml^^' since some of the neural network 
variables are strongly correlated to E^^\ 

b tagging: The cut on the output of the b-tag in the D~ hemisphere is varied 
in the Monte Carlo between 0.43 and 0.57 but kept at 0.5 for the data to 
account for uncertainties in the Monte Carlo modelling of the efficiency of 
the b tagging algorithm. This corresponds to a change in the b tagging 
efficiency of about ±2%. 
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All uncertainties are added in quadrature to determine the total systematic 
uncertainty. 



8 Conclusion 

D*~ — > 7D~ decays are selected in the invariant mass distribution of the photon 
and the D~ meson. The branching ratio of D^" meson decays into ru^ has been 
measured to be 

BR(D7 ^ T-Ur) = (7.0 ± 2.1(stat) ± 2.0(syst)) % , 

in good agreement with the only other direct measurement but with a slightly 
smaller uncertainty [0. From this measurement the D~ decay constant can be 
derived using Eq. [I| and the values in Table |T] to be 

fj3^ = (286 ± 44(stat) ± 41(syst)) MeV 

consistent with theoretical predictions [1-3] for fn^ and with the world average 
fo, = 280 ±48 MeV m. 
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BR(D- 












Gp — 


(1.16639 ±0.00001) X 10"^ GeV~^ 


BR(D*- 


- 7D3-) 


= 0.942 ± 0.025 


\Vcs\ = 


= 0.9891 ±0.016 


BR(r- - 




= 0.1783 ±0.0006 


rrir = 


1.77703 ±0.00003 GeV 


BR(r- - 




) = 0.1737 ±0.0007 


= 


= 1.9686 ± 0.0006 GeV 








^D, = 


(0.496 ±0.01) X 10-12 g 



Table 1: External values used in the calculation of BR(Dg —>■ r u^-) and f^^ [p^ . 
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Source 


A BR/BR CX) 


External Sources: 




/(c - D-) 


20.7 




7.8 


Rc 


2.9 


BR(Dr ^ 7D-) 


2.7 


BR(t H V^Vr) 


0.2 


BR(r~ c'v^Vt) 


0.1 


Monte Carlo: 




background statistics 


11.5 


e(D- ^ r-p,) 


2.3 


e(Ds ^A*) 


0.8 


background rate 


0.3 


background composition 


1.6 


Fragmentation 


1.3 


Lepton spectrum 


0.4 


Detector resolution: 




Tracking resolution 


4.7 


Photon energy 


6.9 


Lepton identification: 




electrons 


2.0 


muons 


4.6 


Neural network inputs 


8.3 


b tagging 


3.0 


total 


28.6 



Table 2: Summary of relative systematic uncertainties on BR(Dg — > r Vr)- 
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Figure 1: Four selected input variables of the neural networks used in the 
electron channel: a) reconstructed energy E-^,^ of the D~; b) output of the neural 
networks trained to find h ^ i decays; c) highest momentum pmax of any particle 
with a charge opposite to that of the lepton in the D~ hemisphere, d) angle 
ttiep.Ds between the lepton in the D~ rest frame and the D~. All distributions are 
normalized to the number of events, Nev The signal contribution in the data is 
about 1% at this stage of the selection. 
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Figure 2: Muon channel: a) output of the neural network against Z — > cc back- 
ground using a logarithmic scale; b) a linear scale, c) output of the neural 
network against Z ^ bb background using a logarithmic scale; d) a linear scale. 
The positions of the cuts are indicated by arrows. 
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Figure 3: Output of the neural network against Z cc background events 
versus the output of the neural network against Z — bb background events. 
The distributions for the muon and the electron channel have been added for a) 
Z — > cc background b) Z — > bb background and c) signal events in the Monte 
Carlo. 
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Figure 4: Distribution of the photon energy for data compared to the Monte 
Carlo simulation. The position of the cut on is indicated by an arrow. 
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Figure 5: Invariant mass m(7D~) of the photon and the D~ candidate for the 
events satisfying all selection criteria. The contributions to the Monte Carlo 
distribution from the signal and from the different sources of background are 
shown separately. The signal region is indicated by an arrow. 
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